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ABSTRACT 

We study accretion induced collapse of magnetized white dwarfs as an origin of 
millisecond pulsars. We apply magnetized accretion disk models to the pre-collapse 
accreting magnetic white dwarfs and calculate the white dwarf spin evolution. If 
the pulsar magnetic field results solely from the flux- frozen fossil white dwarf field, 
a typical millisecond pulsar is born with a field strength ~ 10^^ — lO^^G. The 
uncertainty in the field strength is mainly due to the uncertain physical parameters of 
the magnetized accretion disk models. A simple correlation between the pulsar spin 

and the mag netic field B^, (J^^/IOS-^) ~ (BjlO^^G)-^/^ , is derived for a typical 
accretion rate ~ 5 x W~^MQ/yr. This correlation remains valid for a wide pre-collapse 
physical conditions unless the white dwarf spin and the binary orbit are synchronized 
prior to accretion induced collapse. We critically examine the possibility of spin-orbit 
synchronization in close binary systems. Using idealized homogeneous ellipsoid models, 
we compute the electromagnetic and gravitational wave emission from the millisecond 
pulsars and find that electromagnetic dipole emission remains nearly constant while 
millisecond pulsars may spin up rather than spin down as a result of gravitational 
wave emission. We also derive the physical conditions under which electromagnetic 
emission from millisecond pulsars formed by accretion induced collapse can be a source 
of cosmological gamma-ray bursts. We find that relativistic beaming of gamma-ray 
emission and precession of gamma-ray emitting jets are required unless the dipole 
magnetic field strengths are > lO^^G; such strong dipole fields are in excess of those 
allowed from the accretion induced collapse formation process except in spin-orbit 
synchronization. Strong dipole fields > 10^^ G could in principle be produced in situ. If 
the millisecond pulsars spin up while emitting gravitational wave, the required values 
for the shortest bursts are > lO^^G, pushing the limits of suggested neutron star 
dynamos. 



Subject headings: accretion, accretion disks — binaries: general — magnetic fields — 
pulsars: general — stars: magnetic fields — gamma rays: bursts 
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1. Introduction 

The formation of millisecond pulsars (MSP) is not fully understood. Explosive supernova 
events have been considered as the origin of pulsars (e.g. Mayle & Wilson 1988, Burrows & Hayes 
1996). The high recoil velocities of the newly formed pulsars (due to the off-center explosions), 
however, would likely expel them well away from the Galactic plane in contradiction to the 
observed Galactic MSP distribution. An alternative scheme is accretion induced collapse (AIC) of 
a white dwarf (e.g. Canal & Schatzman 1976, Baron et al. 1987, Nomoto & Kondo 1991, Livio 
&; Truran 1992). However, although AIC is the most probable mechanism for MSP formation 
(Verbund 1993, van den Heuvcl & Bitzaraki 1995), the physics of MSP formation from accreting 
binary systems (Taam & van den Heuvcl 1986, van den Heuvel et al. 1986, Chanmugam & Brecher 
1987, Verbund 1993) and the physical conditions immediately after AIC remain unclear. 

AIC may plausibly lead to pulsars with fast spins in the weakly magnetized white dwarf 
binary systems (Chanmugam & Brecher 1987, van den Heuvel & Bitzaraki 1995), but it is not 
clear how such a mechanism could lead to rapidly rotating and very strongly magnetized neutron 
stars. This is because the strong white dwarf magnetic field is likely to efficiently brake the 
rotation prior to AIC. Narayan and Popham (1989) investigated the consequence of the magnetized 
accretion and the AIC of the white dwarf using the magnetized disk model of Wang (1987). They 
concluded that for an accretion rate M ~ 10^'^ MQ/yr, a MSP would form with a typical limiting 
field strength of ~ 10^^ G. Due to recent improvements in the magnetized accretion disk models 
(Campbell 1992, Cameron & Campbell 1993, Yi 1995, Wang 1995), such qualitative conclusions 
can be improved and further quantified. Some interesting applications of the disk-magnetosphere 
interaction models have recently been made to various magnetized accretion systems (e.g., Yi 
1995, Kenyon et al. 1996, Yi & Kenyon 1996, and references therein), where it is found that there 
could exist a simple spin-magnetic field correlation even when the disk-star system is not in the 
equilibrium spin state (Yi 1995). If such a simple relation exists in the pulsar systems as well, the 
electromagnetic emission from MSPs is likely to reflect such a correlation. 

The spin-field parameter space is important because strongly magnetized MSPs (Bhattacharya 
&L van den Heuvel 1991, Taylor & Stinebring 1986) could be sources of cosmological phenomena such 
as cosmological gamma-ray bursts (GRBs) (e.g. Usov 1992 and references therein) and gravitational 
radiation (e.g. Backer & Hellings 1986, Shapiro k. Teukolsky 1983). Electromagnetically quiet AIC 
could be as frequent as ~ 0.1 — lyr~^ per galaxy which makes the AIC scenario viable for GRB. 
The quiescent nature of the collapse is crucial for this source of GRBs because the bursts have 
not been identified with any precursor activities (for a review, see, Fishman et al. 1994, Fishman 
& Meegan 1995). In the Usov (1992) model, the required conditions for the pre-collapse white 
dwarfs are highly extreme compared with the properties of the known accreting white dwarfs (e.g. 
Wu Sz Wickramasinghe 1993). For a pulsar of surface dipole magnetic field strength ~ 10^^ G, 
the moment of inertia 7* ~ W^^gcm?, the radius i?* ~ 10^ cm, and the spin rate fi* ~ 10^s~^, the 
angular momentum conservation and the flux-freezing lead to a pre-collapse white dwarf with the 
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rotational frequency 

n^^d ~ ^.{h/iwd) ~ io-2s-i (1-1) 

and the magnetic field strength 

B^a - B,{R,/R^d? ~ lO^G (1-2) 

for the white dwarf of radius Ry^d ~ 10^ cm and moment of inertia I^d ~ 10^^(/cm^ (e.g. Frank et 
al. 1992). For the observed cataclysmic variables, such a field strength has not been detected (Wu 
& Wickramasinghe 1993). 

Usov (1992) suggested a possible origin of pulsars with very strong magnetic fields; a binary 
system consisting of two dwarfs. In order for the dwarf binary system to be in mass transfer, 
the binary separation should be close enough to allow the Roche-lobe overflow. Based on this 
observation, Usov (1992) further suggested that the accreting magnetized white dwarf's spin is 
likely to be synchronized with the binary's orbital motion through the direct magnetic coupling 
between the two stars (e.g. Lamb ct al. 1983). Such a scenario could indeed give a pre-coUapse 
white dwarf with ^yjd ^ 10^^. It has been noted, however, that the mass accretion rates required 
for AIC arc as high as > 4 x 10~^ MQ/yr (Livio k Truran 1992, Nomoto & Kondo 1991, and 
references therein). It is not clear whether such high accretion rates could be achieved and 
sustained for an extended period of time in the binary systems where the companion stars are 
dwarfs. It is therefore crucial to examine the pre-collapse conditions of the magnetized binary 
systems. 

If the gamma-ray emission is isotropic, a cosmological GRB (for a review see e.g. Blaes 1994, 
Fishman et al. 1994, Fishman &; Meegan 1995) source must have the luminosity 

« [1.5 X lO^'erg/s] (5 ^ iQ-Ser^/s/cmO {ebkm/s/Mpc) ((l-x)O ^^'^^ 

where is the typical detected gamma-ray flux, Hq is the Hubble constant, x = 1 — (1 + 2;)"^/^, 
2; ~ 1 is the cosmological redshift, and a 1 — 2 is the usual photon index for gamma-ray 
spectra (e.g. Blaes 1994, Fishman Sz Meegan 1995, Fishman et al. 1994, Yi 1993 and references 
therein). It has been noted, primarily based on energetics and time scales, that the MSPs with 
strong magnetic flelds could be interesting cosmological sources of electro-magnetic emission (Usov 
1992, 1994). If the pulsars lose their energy through emission of the gravitational radiation with 
luminosity (Shapiro & Teukolsky 1983), 

= '^me^ . 12 X 10-e.,/.le^ (^)^ (^)' (1-4, 
the characteristic energy loss time scale 

^ - p X i„-3.„- (_^) -\j^y (1-5, 
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could be short enough to account for the observed short durations of the GRBs (Usov 1992, 
Blackman et al. 1996). is the rotational moment of inertia around the 3-axis, fi* is the 
rotational frequency of the pulsar, = 6^2/(2 — 6^2), and 612 is the equatorial eccentricity on the 
plane spanned by the 1-axis and 2-axis perpendicular to the rotation 3-axis, 



ei2 



(1 - {R2/Ri)Y^ (1-6) 



with the radial extent of the pulsar along the three axes taken as i?2, and R3 respectively. For 
the electro-magnetic radiation relevant for GRB, the luminosity (Shapiro & Teukolsky 1983, Usov 
1992,1994, Blackman et al. 1996) 

^™ - ^ - p - io-e../.i {^y {j^y (1-7, 

determines the emission time scale only when ei2 is vanishingly small or Lem > Lgr, where 
-B* = 11^/ is the neutron star's surface magnetic field strength and fi^ is the magnetic moment 
of the star. In order for a pulsar to be a cosmological GRB with short duration and isotropic 
emission, the pulsar has to be both rapidly spinning (fi* ~ 10^s~^) and extremely strongly 
magnetized with ~ lO^^G (Usov 1992). 

However, we find that MSP with > lO^^G magnetic fields, required by the isotropically 
emitting MSP cosmological GRB models, are incompatible with AIC for observed white dwarf 
magnetic fields < 10^ G. If these pulsars do exist, then their fields would have to be generated 
in-situ by a dynamo. Duncan & Thompson (1992) suggest that dynamos can generate lO^^G 
dipole magnetic fields in millisecond or sub-millisecond pulsars because the available differential 
shear and turbulent energy (oc Vil) for conversion to magnetic field is sufficiently large. However, 
complications of imposing such a dynamo are addressed in section 4. Constraints on the MSPs 
formed by AIC and the correlation between pulsar spin and flux frozen fossil field may turn out 
to be most important. Even if dynamos could produce lO^^G fields, we will point out that the 
shortest bursts require field in excess of lO^^G if the emission source is electromagnetic dipole 
radiation and gravitational radiation from MSPs spins up neutron stars. 

Without superstrong fields, cosmological MSP can still be associated with GRB by reducing 
the energy requirements of the field through relativistic beaming. Blackman et al. (1996) 
proposed that a strongly beamed emission from a pulsar with less extreme physical conditions 
could be the origin of the cosmological GRBs. In this scenario, the observed gamma-rays emanate 
from relativistically beamed jets and the observed luminosity could be larger than the intrinsic 
luminosity of the source by a factor due to the relativistic beaming (Yi 1993) where F is the 
relativistic bulk Lorentz factor of a gamma-ray emitting jet (cf. Usov 1994). For a Lorentz factor 
of F > 10^, the required electro-magnetic power (eq. (1-5)) could be lower by a factor of > 10^ 
and hence the required field strength could be < lO^^G for a MSP. In this case, however, the 
characteristic time scale for electro-magnetic emission 



Lem ^ io45^cm2 ) \ l{)^cm ) \ IQi^G ) \ lO^s"! 
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is much longer than the typical GRB durations. Blackman et al. (1996) pointed out that 

precession of the sharply beamed gamma-ray jets is a possible explanation for the observed short 
durations of GRBs. In this picture, the duration of a burst event is determined primarily by the 
time scale on which the jet sweeps by the observer's line of sight. An additional slower precession 
mode, naturally expected in the MSP binary systems (e.g. Thorne et al. 1986), would account for 
the absence of repeaters by keeping the beam from returning to the line of sight. 

When the pulsar rotates with a period P* = 27r/r2* shorter than the critical period (e.g. 
Chandrasekhar 1969) 



e could be as high as ~ 0.1 as a result of the non-axisymmetric gravitational secular instability. 
This result plausibly suggests a source of the non-zero quadrupole moment which is the necessary 
ingredient for the gravitational radiation emission. The short time scale tcR could be the GRB 
duration provided that the pulsar loses its rotational energy, spins down, and Lem decreases as 



emission evolves during the angular momentum loss by the gravitational radiation (e.g. Finn 
&; Shapiro 1990). For the pulsars with fi^ > ^crit-, the application of the conventional dipole 
formula (Ostrikcr &: Gunn 1969, Shapiro & Teukolsky 1983) for the spin-down of the MSPs 
becomes uncertain during the period dominated by the gravitational radiation emission. This is 
because the conventional formula does not take into account the change in the moment of inertia 
(Chandrasekhar 1969). In this sense, it remains to be seen how the MSPs electromagnetic emission 
would evolve right after AIC. Such a question becomes especially interesting for the pulsars with 
> ^crit = '^T^/Pcrit, which is often invoked for the cosmological GRBs (Usov 1992, Duncan & 
Thompson 1992) . 

In this paper, we: (i) Investigate formation of the MSPs from AIC of accreting white dwarf 
binary systems and determine a spin-magnetic field parameter space for the resulting pulsars, 
(ii) Study electro-magnetic emission from rapidly spinning pulsars rotating above the critical 
frequency, (iii) Explore the possibility that even though AIC produced pulsars do not have 
strong magnetic fields, the energy requirement reduction by relativistic beaming can still allow 
AIC pulsars to be associated with cosmological GRBs. In section 2, we describe the magnetized 
accretion disk models. Section 3 gives the results of the AIC and the derived spin-magnetic field 
correlation. In section 4, we examine the electro-magnetic emission from the rapidly rotating 
pulsars and derive the conditions for the gamma-ray emission from the MSPs. We conclude in 
section 5. 




(1-9) 



(Usov 1992, Blackman et al. 1996). It has been, however, unclear how the electro-magnetic 



2. Magnetized Accretion Disk Around a White Dwarf 
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We assume that the white dwarf magnetic field is of the dipole type and is characterized by 
the stellar magnetic flux = B^^Rwd magnetic moment = ^wdRwd where R^d is the 
equatorial radius and Byj^ is the white dwarf surface field strength. The vertical component of the 
magnetic field is approximately given by 

B.{R) = (2-1) 

where we have adopted the usual cylindrical coordinate system {R,(p,z). We consider two 
magnetized accretion disk models in which the accretion disk is threaded by the white dwarf 
magnetic field (Ghosh & Lamb 1979ab, Wang 1987, Campbell 1992, Yi 1995). We briefly describe 
the models with a summary of the relevant formulas. 

The first model assumes that a single turbulent mechanism is responsible for both viscous 
angular momentum transport and magnetic field diffusive loss (Campbell 1992, Yi 1995). In this 

case, using the magnetic diffusivity ijt = avtH where a is the conventional Shakura-Sunyaev 
viscosity parameter (Frank et al. 1992) and H is the vertical disk scale height, the azimuthal 
component of the field is given by 

B^R) = l^halzIlliB.iR) (2-2) 
a I Ik 

where Qyj^ is the white dwarf rotational frequency, ^Ik = {GM,^d/ R'^Y^'^ is the Keplerian disk 
rotational frequency and 7 ^ \{R/n){dVL/ dz)\ measures the vertical velocity shear length scale 
between the Keplerian disk midplane and the corotating stellar magnetosphere (Wang 1987, 
Campbell 1992, Yi 1995). 

The larger the ratio 7/a, the larger the diffusion of the external stellar field into the disk 
becomes (Ghosh & Lamb 1979ab, Yi 1995). The torque on the star is given as 

v_7 ,, l-(8/7)(i?o/i?e)^/^ 

where A^o = M{GMyjdRoY''^ and i?c = (GM^„rfP2^/47r2)V3 ig the corotation radius with 
Py^d = '^T^/^wd (Yi 1995, Wang 1995). The disk inner radius Rq below which the disk is 
magnetically disrupted is determined by 



RcJ 



1 



RoV^' 
R, 



(2-4) 



where A = 2{-i/a)BlRl/ M{GM^dRcY''^ and = Hwd/Rl (Wang 1987, Yi 1995). Eq. (2-3) is 
the combination of the three physically well-defined torque components: In the region between Rc 
and i?o, the disk's Keplerian rotation is faster than the stellar rotation and the disk exerts the 
spin- up torque on the star. In the region outside i?cj the disk rotation is slower than the stellar 
rotation, which results in the spin-down torque. The disk material accreted by the star adds the 
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specific angular momentum (GMyj^Ro)^^^ ■ Using the characteristic parameters, A becomes 



(Yi 1995). The torque vanishes when Rq/Rc = 0.9148. As Rq/Rc —>■ 0, the torque asymptoticahy 

approaches A'^ 7Nq/Q. 

In the second model (Aly & Kuijpers 1990, Livio & Pringle 1992, Wang 1995), it is explicitly 
assumed that the magnetic field pitch (the ratio of toroidal to vertical field strength = \B(f,/Bz\) 
is limited by the wind-up of the stellar field lines due to the velocity shear between the Keplerian 
disk and the stellar magnetosphere corotating with star (cf. eq. (2-2)). This could be the case if 
the magnetosphere is nearly force-free (i.e. V x B oc B) and the winding-up of the magnetic field 
lines is limited by reconnection (Aly & Kuijpers 1990) within the stellar and disk magnetosphere. 
The details of the magnetosphere arc hard to model. But, if the reconnection occurs on a time 
scale of order of the vertical shear time scale between the star and the disk (cf. 7), then the pitch 
of the field will be limited to a factor of order unity (Aly & Kuijpers 1990, Livio k. Pringle 1992, 
Wang 1995). In this case the azimuthal field is 

B^{R) = 7„„,^^^S,(i?) (2-6) 

for R < Rc 01 iluid ^ (Wang 1995). The parameter jmax sets the maximum pitch of the field 
configuration. For R> Rc or Clwd > ^K, we get 

B^R) = lmoJ^^^^^^B,{R). (2-7) 

We note that ^max of this second model replaces the ratio 7/a in the first model for R > Rc- 
The quantity ^max mainly affects the region outside R^ where in the first model the magnetic 
pitch increases as ~ VL^fi/Q,K oc R^/"^. Practically, we will see that the difference between the two 
models is not serious because the magnetic torque in the outer region rapidly decreases due to the 
steep decrease of the dipole field B^ oc R~^. 

The inner truncation radius Rq is determined by the same equation as for the first model, 
with 7/0; replaced by ^max in the expression for A. We have 

MiGM^aRcY/^ 

in eqs. (2-4), (2-5). The expression for the torque on the star becomes (Wang 1995) 

7 1 - {8/7){Ro/R,y/^ + {2/21){Ro/Rcr 

which vanishes {N = 0) when Rq/Rc = 0.9502. 
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The difference between tfie two models is sfiown in Fig. 1 using tfie torque behavior as a 
function of the disk truncation radius scaled by the corotation radius. Note that, as mentioned 
above, the difference between the two models is not significant because of the rapid decrease of 
the dipole field for R > R^.. We therefore expect small differences in the spin evolution of the 
white dwarf provided that the binary system parameters are similar. We take a = 0.1, 7 = 1, and 
Imax = 1 and discuss the uncertainties related to these chosen parameters later. 

When the mass accretion rates are large and the magnetic fields are weak, Rq becomes 
comparable to the stellar radius Ry,d and the accretion disk extends down to the stellar surface. 
In this limit, we use the purely hydrodynamic accretion torque M {GMyjdRwd)^^"^ in place of the 
magnetic accretion torque (eqs. (2-3), (2-9)). Note that in non-magnetized accreting systems, 
sustained accretion is possible even with spin-down when the star spins near the break-up rate 
~ (GM^dZ-R^d)^^^ (Popham & Narayan 1991, Paczynski 1991). We refer to objects that reach 
this limit as critical rotators (Narayan Sz Popham 1989). We, however, do not consider these cases 
in detail because their direct collapse to strongly magnetized neutron stars is not likely. 



3. Non-Explosive Accretion Induced Collapse and Neutron Star Spin 

3.1. Model Parameters for Accretion Induced Collapse 

Initial Mass of White Dwarf: We assume that the collapse of an 0-Ne-Mg white dwarf occurs 
when the central density reaches the electron capture threshold for Mg at pc = 3.16 x IQ'^g/cm^ 
(cf. Nomoto & Kondo 1991, Shapiro & Teukolsky 1983). The critical Chandrasekhar mass of the 
white dwarf is taken as M^^ = 1.39M0. 0-Ne-Mg white dwarfs of mass My,d = 1.2 - 1.37M0 

can form from 8 — IOMq progenitor stars (Nomoto & Hashimoto 1988). We therefore choose 
M.fjj(i = 1.25M0 and M-^jd = 1.35Af0 as initial masses of the white dwarfs (see also Livio & Truran 
1992). It turns out that within this range the initial mass only significantly affects the time scales 
to reach collapse for a given accretion rate but not the field-spin relation (sec below). 

Mass accretion Rate: In principle, the nature of the collapse depends on the accretion rate, 
but the range of accretion rates, M, for AIC has been constrained. Low mass accretion rates, 
M < 10~^ MQ/yr, are believed to lead to novae or similar eruptions (Paczynski and Zytkow 
1978), which would decrease the mass of the accreting white dwarf (Livio &; Truran 1992). 
Intermediate mass accretion rates, lO^^MQ/yr < M < 4 x 10~^MQ/yr, likely lead to off-center 
He detonation (Nomoto 1987, Nomoto & Kondo 1991). This leaves high mass accretion rates 
M > 4 X 10~^MQ/yr = 2.5 x lO^^g/s as the most likely window for AIC (Livio & Truran 1992). 

For the plausible mass range of the 0-Ne-Mg white dwarfs and an accretion rate, 
M ~ 10~^ MQ/yr of observed white dwarf binary systems, the accretion rate has to be sustained 



-9- 



for > lO^yr. Such accretion rates are generally expected from a mass-losing lower giant branch star 
of mass ~ IMq that climbs up the giant branch (Webbink et al. 1983). It is questionable whether 
such a high accretion rate could be achieved for a dwarf binary system where the secondary star 
in the binary has a mass <C IMq (cf. Usov 1992). For most of the observed white dwarf binary 
systems (cataclysmic variables), the mass accretion rates generally appear to be much lower than 
the typical Eddington rate M < 1.4 x lO^^g/s (e.g. Prank et al. 1992). This may imply a narrow 
range of allowed mass accretion rates for pre-AIC binary systems. Although very high mass 
accretion rates may also lead to collapse, they would be very rare and the duration of such activity 
in the binary systems is most likely to be short (e.g. Smak 1984). Wc therefore mainly consider 
mass accretion rates near M = 5 x 10~^ MQ/yr. Although we consider a larger rate {W~'^ Mq/uv) 
and a smaller (2 x 10~^ MQ/yr) rate for comparison, the differences do not affect our conclusions. 

Rotating White Dwarfs: Our calculations require a rotational law, radial density profile, 
and moment of inertia. Most importantly, we need a mass-radius relation as the mass increases 
through accretion and the star spins up. We specifically adopt the results in Hachisu (1986) which 
are relevant for uniformly rotating white dwarfs (e.g. Narayan &; Popham 1989). Here the white 
dwarf as a whole is approximated as a rigidly rotating solid body. When the white dwarf is rapidly 
spinning, the eccentricity of the star with respect to the rotation axis becomes significant. The 
results of Hachisu (1986) can be interpreted essentially as a family of solutions corresponding to a 
combination of parameters, Myj^, angular momentum Ju,d, Rwd, and R3/R1 (cf. eq. (1-6)) for a 
given central density of the axisymmetric star. 

Spin Evolution: When the white dwarf mass increases significantly through accretion and 
the radius of the rapidly rotating star varies in response, the spin of the star is affected both 
by the accretion disk torque and by the change in the moment of inertia. Angular momentum 
conservation gives 

-IT = ^""'^ + ^""^^ = ^ ^^-^^ 

where J^^ = Iwd^wd = P^wdRwdf^wdj the moment of inertia I^id = ^M^^^^'^d ^^'^ constant 
P = 0.08 depends on the details of the white dwarf structure including the rotational fiattening 
and the radial density profile (Hachisu 1986). Using this expression, we get 



dt 
or 



N 1 dM^d _ 2 1 dR^d 



(iMwdRwd^wd Myjd dt Ryjd dt 



d^wd Q 

~dr ~ 



1 dMy^d _ 2 1 dRy,d 

Jwd Mwd dt Rwd dt 



(3-2) 
(3-3) 



We further assume that the magnetic flux threading the stellar surface is conserved, i.e. 

Bwd = Bwd,i f (3-4) 



■wd 



where the subscript i denotes the initial epoch at which the accretion and spin evolution begin. 
Pinally, we assume that all white dwarfs spin sufficiently slowly at the beginning of their evolution. 
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so that, due to the short spin-up (-down) time scales, the details of the initial spin rates are 

not necessary (e.g. Yi 1995). Once the constant mass accretion rate is chosen and the initial 
field strength is assigned, the mass increase and torque arc calculated while the radius and the 
magnetic field of the white dwarf arc updated accordingly. 

There are three possible cases for collapsing white dwarfs (e.g. Narayan and Popham 1989): 

(i) When the star rotates faster than the critical rate, roughly equal to the equatorial 
Kcplcrian rotation ~ (GM^^/i?^^)"*^/^, the star can lose its mass due to the centrifugal force (i.e. 
equatorial mass shedding) and the rotational support prevents any further collapse. If the mass is 
continuously added to the star with the critical rotation, it is unclear how the star would respond. 
In the non-magnetic, steady state case, the purely hydrodynamic accretion can proceed with 
continuous spin-down as shown by Popham and Narayan (1991) and Paczynski (1991). In the 
magnetized accretion case, the accretion occurs mainly along the polar accretion column (Frank et 
al. 1992). This polar accretion may continue despite the loss of matter near the equatorial plane. 
Given this complex situation, we simply assume that the equatorial mass loss and the polar inflow 
balance, so that the net accretion essentially ceases when the star reaches the critical rotation. 

For comparison with other possible evolutions, we calculate the hypothetical MSP spin and 
magnetic field after a flux freezing and angular momentum conserving collapse of the white dwarf 
at the critical rotation. Practically, these objects correspond to those with weak magnetic fields, 
which are about ~ 80% of the known cataclysmic variables (Morris et al. 1987), and hence they 
are unlikely to be a source of intense electromagnetic emission. Such MSPs, if they indeed form 
from AIC, could be a significant source of the gravitational radiation. 

(ii) Consider now the case for which the white dwarf reaches the collapse density and the 
angular momentum of the star is below the critical rotation limit but above ~ 6 x W^^gcrv? /s. 
Here the collapse of the star is likely to be halted after shrinking to a size for which the centrifugal 
forces prevent further collapse. These objects are often referred as fizzlers (Shapiro and Lightman 
1976, Tohline 1984). The critical angular momentum for the fizzlers is taken approximately as 
J* = ~ 6.3 X IQ^^gcm? /s for P* = 27r/J7* 10~^s. Fizzlers with fast spins are likely to 
be unstable to non-axisymmetric perturbations (Lai &: Shapiro 1995). The resulting triaxial 
configurations could lose angular momentum through gravitational radiation which could push the 
object to below the critical angular momentum limit. The collapse would then resume leading to 
the formation of a pulsar. Below, we treat the collapse of fizzlers as direct collapses by neglecting 
possible angular momentum loss at the fizzler stage (cf. Durisen et al. 1986, Williams & Tohlin 
1988, Houser et al. 1994). In this case, we also expect relatively small magnetic fields for most of 
the fizzlers and these cases have only limited relevance for strong electromagnetic emission and 
any associated cosmological GRBs. 

(iii) Finally, when the angular momentum is below the fizzler limit and the central density 
exceeds the critical value pc ^ 10^'^gi/cm^ (Nomoto &: Kondo 1991, Shapiro & Teukolsky 1983), 
the stars could collapse directly to neutron stars. In this case, the field strength is expected to be 
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large and the rotation is relatively slow. We assess how likely it is to get rapidly spinning strongly 
magnetized neutron stars in this case. We calculate the final spin and magnetic field assuming 
conservation of angular momentum and magnetic flux. The magnetic flux of the collapsed object 
is then given by i?* = {Ru,d/R*)^l^wd and the rotational velocity = {I-wd/ I*)^wd where we take 
= 10^ cm and /* = 10^^ gcm^. 



3.2. Correlation between Spin and Magnetic Field 

Fig. 2 shows the main results of the spin evolution. In all panels, the dashed lines correspond 
to the cases where the white dwarfs reach the critical rotation before collapse. The field strengths 
and spins for the critical rotators shown in Fig. 2 are hypothetical values (strictly for comparison) 
based on the assumption that the stars begin to collapse at the critical rotation. These objects are 
characterized by fast spins and weak magnetic fields reflecting the inefficient magnetic braking of 
the disk-star interaction. The dotted lines correspond to the fizzlers which are assumed to collapse 
down to the neutron stars after reaching the critical central density. Since the fizzlers arrive at the 
collapse line with more angular momentum than the objects which collapse directly (solid line), 
the neutron stars formed from fizzlers would be observed as those with somewhat higher spins 
and weaker magnetic fields. In our models, the resulting magnetic fields are not strong enough for 
the fizzlers to be strong emitters of electromagnetic radiation but they could be strong sources of 
gravitational radiation. 

In all panels, thin lines are the results of the pitch-limited magnetized disk model (eq. (2-9)) 
and the thick lines are for the model with the diffusive magnetic field loss (eq. (2-3)). In the 
upper panels, the initial mass of the white dwarf is M^d = 1.25Mq. ^From (a) to (c), the mass 
accretion rate increases from 2 x 10~^ MQ/yr to W~'^ Mq/ut. The highest accretion rate results 
in the fastest spin. We note that the highest mass accretion rate used in (c) has not been seen in 
typical white dwarf binary systems such as the cataclysmic variables (Frank et al. 1992) except 
during brief outbursts (Smak 1984). Although some transient phases of high mass accretion rates 
cannot be ruled out, such high accretion rates are less likely to be sustained for an extended period 
of time. 

In the lower panels (d)-(f), the initial mass of the white dwarf is increased to M^d = I-SSMq. 
The differences between the upper and lower panels, due to the difference in the white dwarf 
mass, are significant only for the critical rotators. We observe that the typical high spin pulsars 

~ 10^s~^ would be born with 5* ~ lO^-'^G. This result is qualitatively consistent with the 
result of Narayan &; Popham (1989). The derived correlation between spin and magnetic field 
strength is one of the main quantifiable characteristics of the MSPs formed from AIC. Different 
choices of the magnetic field model and the initial white dwarf mass only slightly affect the results. 

Since there is little difference between results in the two different white dwarf initial mass 



- 12 - 



cases, we can write down a simple spin-magnetic field correlation for a given mass accretion rate. 
The two distinctive magnetic field models wc discussed in section 2 also do not result in serious 
numerical differences. For the diffusive loss case (thick lines) and for the mass accretion rate 
M = 5 X 10~^ Mq/ut, the results in Fig. 2 show that the spin-field correlation gradually varies 
from cx B~^'^^ near the fizzler regime to f2* oc toward the slow spin regime. The latter 

correlation is precisely what is expected in the spin equilibrium situation (cf. eqs. (2-3), (2-5)) 
before AIC; the equilibrium relation is self-similar because both the field and the angular velocity 
vary inversely with B?. 

The results indicate that the spin equilibrium is not achieved for most of the rapidly spinning 
white dwarfs. This is not surprising given the fact that the mass and moment of inertia are 
constantly changing during the pre-collapse accretion phase. For MSP (or sub-MSPs, i.e. near 
(i.e. near ~ 10^s~^), and moderately strong magnetic fields, the spin-field correlation is best 
described by 

where we again note that the two initial masses have little effects on the cases of the fizzlers and 
direct collapses. The initial white dwarf mass only affects the early spin evolution to critical 
rotation before collapse. 

For the pitch-limited torque model (thin lines) and for the same parameters as above we get 

This simple spin-field correlation is one of the main results of the paper. An analogous spin-field 
correlation has also been predicted in the magnetized accretion model for the protostellar systems 
(Yi 1995, Kenyon et al. 1996). 

There are two reasons for the small differences between the two torque model predictions. 
The first is that in model two, the assumed maximum pitch -jmax effectively reduces the spin-down 
torque contributed by the region outside the corotation radius. This effect is less important as the 
dipole field strength rapidly decreases outside the corotation radius and the torque contribution 
from the outside region is relatively small. The second reason is due to the scaling in the effective 
measure of the magnetic field strength. For the diffusive loss model, we have adopted a = 0.1 and 
7=1 which combines to give 7/a = 10. As we have mentioned earlier, ^max essentially has the 
same physical effect as the combination 7/a. Therefore, for our choice of 7/a = 10, the effective 
field strength in the diffusive loss model is roughly three times larger than that of the magnetic 
pitch-limited model. The results in Fig. 2 reflect this effect and show smaller spin rates for the 
same stellar field strength and the mass accretion rate. 

Although the spin equilibrium {N = 0) condition is not always reached during the pre- AIC 
spin evolution shown in Fig. 2, the equilibrium spin provides a useful scaling for the derived 
spin-field correlation and especially for the dependence on the model parameters 7/a or 7^00; • 
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Assuming that the final pre-coUapse white dwarf spin period is the equihbrium spin period, using 
eq. (2-5), we can show that the pre-collapse white dwarf parameters have to satisfy 

a)\WG) \W^) [m^] llOOij [wWs) ' ^'"'^ 

which shows the equihbrium spin-field correlation Pyjd oc B^J as noted earlier. For the 
representative parameters Ryj^ ^ 10^ cm, My,^ ^ 1.3M(r), ^ 10~^s~^ relevant for AIC, 



1/2 



which is consistent with the derived result, eq. (3-5), within a factor of order unity. That 
is, assuming fiux-freezing and angular momentum conservation for AIC, the derived spin-field 
correlation eq. (3-5) can be used in eq. (3-7), which results in a relation very close to eq. 
(3-8). Taking into account the anticipated uncertainty in ■y/a, for instance, by about an order of 
magnitude, our derived field strengths would be uncertain at most by a factor of 3. Typically, 
for the cataclysmic variables a ~ 0.1 is often quoted (Frank et al. 1992). 7 could in principle 
range from ~ 1 to 



where H is the vertical thickness of the accretion disk (e.g. Prank et al. 1992). This indicates 
7 ~ 0(1) for high mass accretion rates M > 3 x W^^g/s appropriate for AIC. In the case of the 
pitch-limited case, from eq. (2-8), we arrive at a result similar to eq. (3-8) with "jmax in place of 
7/0;, i.e. 

which is also consistent with eq. (3-6) within a factor of order unity. We note that "fmax is 
expected to remain close to unity as long as the reconnection time scale is close to the time scale 
of the vertical velocity shear (Aly &; Kuijpers 1990). Eqs. (3-8) and (3-10) are consistent with 
the dependence of the spin-field correlation on 7/0; or "fmax seen in eqs. (3-5) and (3-6). We 
conclude that the numerical difference between the two spin-field correlations is largely due to the 
different effective strength of the magnetic field determined by 7/a or jmax- It is gratifying that 
the uncertainty in 7 and 'jmax is likely to be at most at the level of order unity. 



3.3. Star-Star Direct Magnetic Coupling 
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We now briefly examine the possibility of magnetic coupling between the accreting white 

dwarf and the mass-losing secondary. Such magnetic coupling between pre-collapse, close binary 
stars has to be considered because it may lead to a synchronization of the white dwarf's spin with 
the orbital motion (Lamb et al. 1983, Lamb and Melia 1987). As pointed out by Usov (1992), the 
short spin period of the pre-collapse white dwarf could be due to this type of magnetic coupling. 
A critical assessment in the context of AIC is necessary. 

The magnetohydrodynamic torque between the white dwarf and the orbiting secondary is 
estimated to be (Lamb et al. 1983) 

Nun ~ (t>lbDRl{^^^d|D^f (3-11) 

where (/> < 1 is the fractional area of the secondary star threaded by the magnetic flux, 75 is the 
pitch of the magnetic fields connecting the two stars, D is the binary separation, and R2 is the 
secondary star's radius. Assuming 7b ~ 1 (Low 1982, Aly 1984), we can estimate the importance 
of the magnetic torque due to the secondary. Once the system is synchronized, we expect 7;, = 0. 
If the secondary star is magnetically active, ^ ~ 1 is expected and we can use this as an upper 
limit. We estimate 

Nun ~ [3 X lO^^G^cm^^l (07^) ( ~' ( \—] ' ( J!^"" A ' (3-12) 

where jiwd — BwdRid ~ W'^^Gcm^ is the magnetic moment of the white dwarf with R^^ 10 cm 
and Byjd ~ lO^G. We note however that /ly^d could be as high as lO^^Gcm^ for very strongly 
magnetized white dwarfs with Byj^ ~ lO^G as suggested by Usov (1992). Although such white 
dwarf magnetic fields ~ lO^G have not been directly observed, the possibility that a small fraction 
of magnetized white dwarfs may have such strong fields cannot be ruled out. We compare N^in 
with the typical magnetic torque due to the accretion disk 



AT ~ [1 X 1036G2^m3] ( — ^ 11^1 ( , 7" ) (3-13) 

(cf. eqs. (2-3), (2-9)). A/ftj„ becomes comparable to N only when D <^ lO^^cm or Byj^ S> lO^G. 

In order for the accretion through the Roche lobe flow to occur, the binary separation D has 
to satisfy 

5 \l/3 



D<m(jf-J (3-14) 

where M2 is the mass of the secondary star and the binary mass ratio q = M2/Myj(i (Frank et 
al. 1992, Paczynski 1971). For N^n > A^ or D< lO^^crn, M2 <C 0.7[g/(l + q)]'^-^^ . We have 
assumed the dwarf mass-radius relation for the secondary together with the binary orbit formula 
originally due to Paczynski (e.g. Paczynski 1971). We do not consider massive companions due 
to their short evolution time scales. The synchronization of the white dwarf spin with the binary 
orbit could occur on a time scale 

tsync Iwd^wd/^bin 
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(3-15) 

The white dwaxf spin would be synchronized with the binary orbital frequency 



1.5M 







5 X lO^^cm 



(3-16) 



which is interesting only when the secondary star is a rather extreme dwarf star. Otherwise, the 
accretion and collapse would have taken place by the time the synchronization could have taken 
place. Assuming that the accreted mass is typically ~ O.IMq for collapse, which we take as the 
minimum mass for the secondary, the typical size of such a dwarf star would be R2 ~ lO^'^ cm 
using the dwarf mass-radius relation. That is, unless the binary stars are in physical contact (i.e. 
a contact binary) , the torque from the accretion disk dominates in white dwarf binary systems 
with less extreme conditions. We conclude that the accretion disk torque is dominant except when 
the binary separation D < lO^'^cm for M > 10^^ g/s provided that white dwarfs have magnetic 
fields not much stronger than the typical field strength ~ a few xlO^G. 

If the binary stars get close enough to have a short synchronization time scale, tgync, the 
orbital evolution due to the gravitational radiation emission is not negligible. The orbital evolution 
time scale due to the gravitational radiation emission is 



^orb,GR 



1 dPorb 



orb 



dt 



which is shorter than the synchronization time scale unless the binary separation D is much 
shorter than ~ 5 x 10^° cm. For B^d ~ lO^G, the synchronization due to the direct star-star 
magnetic coupling is most likely to operate in systems with M <^ W^^g/s, which is not relevant 
for AIC. We note however that the condition tsync < torb,GR could be easily met if white dwarfs 
are very strongly magnetized with Byj^ ^ lO^G (eqs. 3-15,3-17). In such a case, N^in > N (eqs. 
3-12, 3-13) and the spin-orbit synchronization is achieved. 

Observationally, the spin-orbit synchronization occurs in AM Her systems, which are a 
subclass of the cataclysmic variables with relatively strong magnetic fields ~ lO^G. The typical 
white dwarf spin periods ~ lO^s would certainly be too long to produce P* ~ 10~^s in the event 
of AIC. The spin periods shorter than ~ lO^s are mostly observed in the DQ Her systems which 
are relatively weakly magnetized with B^^ < lO^G, which indicates that the short spin periods 
arc more likely in systems with accretion disks and weak fields (Warner & Livio 1987, Warner 
& Wickramasinghe 1991). This observation is consistent with our estimates on the relative 
importance of the accretion torque, eqs. (3-12), (3-13). As noted by Chanmugam &; Brecher 
(1987), the observed cataclysmic variables have no trouble being candidates for the MSPs through 
AIC. In this sense, the role of the star-star synchronization is not essential in producing MSPs. 
The known number of cataclysmic variables ~ 10^ in our Galaxy (de Cool & van Paradijs 1987) 
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indicates that only a small fraction of them are required to undergo AIC to produce ~ 10^ low 
mass X-ray binaries with neutron stars (Chanmugam &; Brecher 1987). 



4. Electromagnetic Emission from Rapidly Rotating Pulsars 

4.1. Electromagnetic Emission and Relativistic Beaming 

Our derived correlation between the spin and the magnetic field strength indicates that the 
electromagnetic dipole emission from MSP immediately after AIC has a very simple scaling. The 
dipole emission formula eq. (1-5), 

Lem oc BlQ.i (4-1) 

leads to a power-spin correlation 

LEMO^ B^nt oc 57^^^ oc nl^^ (4-2) 

where we have used the spin- field correlation cqs. (3-5), (3-6). Taking into account the numerical 
factors in eqs. (3-5) and (3-6), we write the single scaling for the spin-field correlation 



lO^s-i; V7xl0iiG 



where -Be// = <5^* with 5 = \fy/a or 5 = ■s/'jmax depending on the model. The electromagnetic 
power from the pulsars is then given by 



LEM-ilO'^ergMSl-^] ( ) . (4-4) 



W^cmJ VlO^ 



This is much smaller than that required in eq. (1-3) indicating that GRB emission would require a 
high degree of beaming or anisotropy (Blackman et al. 1996); otherwise the burst event would not 
be visible at cosmological distances. Assuming a nominal gamma-ray emission efficiency ^ = 0.1, 
the observed gamma-ray luminosity L-y > 10^^ erg/ s would require a degree of the relativistic 
beaming characterized by the bulk Lorentz factor T, 

r = {L^/^Lem)'/'' > 104<5V2(^/0.1)-V2 (4.5) 

for a typical MSP with $7* ~ 10^s~^ and i?* ~ 10^ cm. This is the first necessary ingredient of the 
present cosmological GRB scenario. 

The low intrinsic power also presents another problem as pointed out in Blackman et al. 
(1996). Using the expected electromagnetic power, eq. (4-4), the characteristic duration of the 
electromagnetic emission is (cf. eq. (1-8)) 

^ ~ f5 X 10^.1 (^Y' ( J* A (^V^' 



= ~ [5 X 10«s] TT^if^ (4-6) 
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if the pulsar's rotational energy is lost only through the dipole radiation. When the pulsar's 

rotation is fast enough to trigger the gravitational secular instability (Chandrasekhar 1970ab), the 
pulsar's quadrupole moment can be significant and the pulsar's angular momentum will be lost on 
a time scale ~ tcR (eq. (1-5)). If /* remains constant during the emission of the gravitational 
radiation, the characteristic time scale of spin-down could be ~ toR on which Lem oc il^ would 
also rapidly decrease. This time scale, eq. (1-5), which depends sensitively on e ~ 0.1, is short 
enough for the short duration cosmological GRBs (Usov 1992) if rapid spin-down occurs as a 
consequence of gravitational radiation emission. The change of inertia as a function of angular 
momentum is however far more complex than the simple case of = constant (e.g. Chanrasekhar 
1969, Finn & Shapiro 1990) and the spin evolution during the gravitationally unstable phase may 
cause spin-up rather than spin-down. In this sense, it is not clear how Lem would evolve and 
what the characteristic time scale for the evolution of Lem is. We therefore must ask: what is the 
characteristic time scale for variation of the electromagnetic power during the evolution driven by 
gravitational radiation emission? 



4.2. Magnetized Neutron Stars as Homogeneous Ellipsoids 

We address this question by approximating the pulsars as magnetized homogeneous ellipsoids 
(Chandrasekhar 1969). When a homogeneous ellipsoid rotates with a period shorter than the 
critical period eq. (1-9), the axisymmetric Maclaurin spheroids bifurcate into two branches (e.g. 
Chandrasekhar 1969, Shapiro & Tcukolsky 1983). In addition to the continuous sequence of 
the Maclaurin spheroids, non-axisymmetric ellipsoids also exist as an equilibrium configuration. 
Maclaurin spheroids with faster than critical spins are unstable to non-axisymmetric secular 
instability and triaxial Jacobi ellipsoids have often been used to model the equilibrium configuration 
resulting from this instability. The non-zero quadrupole moment of a non-axisymmetric star is 
the origin of the gravitational radiation. As shown by Chandrasekhar (1970ab), the equilibrium 
Jacobi ellipsoids can spin up even though angular momentum is radiated away. This will strongly 
affect the electromagnetic emission because the dipole emission formula involves the magnetic 
moment and the spin rate, not the total angular momentum. When the pulsar rotates at a rate 
below the critical rate, the non-axisymmetry is irrelevant and the electromagnetic emission gives 
a long time scale (eq. (4-6)). But, if the pulsar rotates faster than critical, the pulsar's overall 
electromagnetic emission time scale has to be separately considered. 

A large non-axisymmetry can also result from deformation by strong enough magnetic fields 
(Ostrikcr &: Gunn 1969, Usov 1992). This may affect the properties of pulsars' equilibrium 
configuration and may change the bifurcation point at which the Jacobi ellipsoids branch out and 
hence the eventual gravitational radiation emission (Tsvetkov 1983, Bonazzola & Gourgoulhon 
1996). When the magnetic fields are very strong, the equatorial eccentricity caused by the 
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asymmetric magnetic tension is estimated as 

where Bp is the poloidal field component and iJ^ is the toroidal field component. Therefore, for 
the fields we have derived ~ 10^^~^^G, such deformation effects should be negligible. The pulsar is 
treated as a uniformly rotating homogeneous ellipsoid. We consider the dynamical evolution of the 
ellipsoid and assume that magnetic field responds passively to the change in pulsar configuration. 

Following Chandrasekhar and Esposito (1970), the energy loss rate of the pulsar due to 
gravitational radiation is given by 

dE^^__G_/d^\/^\ 

dt 45c5 \ dt^ / \ dt^ I ^ ' 

and the angular momentum loss is given by 

dJ, AG /d^I,kd^Iik\ 



dt 5c5 \ dt^ dt^ 



(4-9) 



where E^, is the total energy of the pulsar, J* is the angular momentum, Dij = 3/jj — 6ijlkk, and 
lij is the moment of inertia tensor given by the standard formula (Chandrasekhar 1969) 

lij = / p^XiXjd^'x (4-10) 

Jv 

where is the uniform stellar density and x is the displacement vector from the center. We have 
adopted the standard tensor notation with i,j,k denoting coordinate indices. The () denotes 
averaging. 

For a quasi-steady pulsar in the frame rotating with angular velocity fi*, eqs. (4-8), (4-9) 
become 

dE^ _ ^ ,2 



dt 5c5 



dt 5c5 

which recovers the familiar equation 



{hi - h2Y (4-11) 



{hi - h2? (4-12) 



1 dE^ dJ^ 



dt dt 

The triaxial Jacobi spheroids have axes i?2, and i?3, related by the constraint equations 

(R1R2V ^ (A3_ 
\ Rs J U12 



and 

■>2 



(4-14) 

ni = 25i2 (4-15) 
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where 



Ai2=RlR2R3j^ (^2+^)3/2(^2 + ^)3/2(^2+^)1/2 

dr 



— R1R2R3 



B12 — R1R2R3 



POO 

Jo Wl 



2 + r)V2(i?2+^)l/2(^2+^)3/2 

rclr 



(4-16) 
(4-17) 
(4-18) 



/o (i?f + r)3/2(i?2+^)3/2(^2+^)l/2 

and we have taken the 3-axis as the rotation axis. Another constraint comes from the assumed 
homogeneity and incompressibihty, i.e. 



R1R2R3 = R = constant 



(4-19) 



We define the axis ratios I2 = R2/R1 and = R-^/Ri from which we get Ri = Rihh) The 
coupled differential equations for the evolution become 



C3C5 



dr CiCry — C2C4 



dh 



C3C4 



dr C2C4 — C1C5 

where r = t/tc is the dimensionless time in units of the characteristic gravitational time 
(Chandrasekhar 1970ab) 



(4-20) 
(4-21) 



to 



25 I R_\ R 
18 i^J 7 



[6.6 X 10-^ s] 



R 



I AM, 



0, 



106 



cm 



(4-22) 



where Rg = 2GM^/(? is the Schwarzschild radius. The relevant quantities used in eqs. (4-20), (4-21) 
are 



-(llli-l)Xi- 3(1 + ^2)^2 



C2 = -l2{l + ll) (^^1+^: 

C3 = -{hl3f'\l-ll?Xl 



Ca — I2 \ 3^2 — InX- 



2^3 



/2 

'2 



with 



Xi — I2I3 



X2 — hh 



-I3 ^3X5 — 

rdr 



(l+r)3/2(Z2 + r)3/2(;2+^)l/2 

rdr 



(l+r)3/2(Z2 + r)5/2(;2+^)l/2 



(4-23) 

(4-24) 
(4-25) 
(4-26) 

(4-27) 

(4-28) 
(4-29) 
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/"OO vdv 

/■OO vdv 
^'"^'^'io (l + r)V2(/2+^)i/2(;2+^)5/2- (4-32) 

For the numerical integration, as in Chandrasekhar (1970ab), we adopt an initial equilibrium 
configuration given by h = R2/R1 = 0.43223, I3 = R3/R1 = 0.34506, 2Xi = ^l/irGp^ = 0.28403, 
which corresponds to the stable equilibrium configuration on the sequence of the Jacobi ellipsoids. 
^From the time evolution, the stellar rotation can be read using f2*/(7rGp*)^/2 = (2X1)^/^. 

The result of the integration, shown in Fig. 3, is similar to that of Chandrasekhar (1970ab) 
except for the overall evolution time scale. Wc summarize the main result as follows: (i) The 
rotation of the pulsar changes from a;* = 0.5329 (initial) to lo^ = 0.612 (final) at which point 
the pulsar has transformed into the Maclaurin spheroids (Fig. 3(a)). (ii) The equatorial area 
of the pulsar is 'KR1R2 = 7rR'^{l2/ll)^^^ ■ The ratio between the initial and final equatorial areas 
is 1.537/1.433 = 1.073 reflecting the deformation due to the change of pulsar rotation, (iii) The 
3-axis eccentricity, 613 = (1 - ll)^^^, changes from 0.9386 (initial) to 0.8127 (final), (iv) The 
equatorial eccentricity, ei2 = (1 — ^2)^^^! decreases from 0.9018 toward e = continuously (Fig. 
3(c)). (v) The gravitational radiation luminosity decreases roughly on a time scale O(100)tG as 
expected (Fig. 3(d)). The evolution of the Jacobi ellipsoid, due to the emission of the gravitational 
radiation, spins- up the pulsar as was noted by Chandrasekhar (1970ab). The spin-up time scale is 
~ O{10^)tc which is similar to tcR in eq. (1-5). 

Despite the spin-up of the pulsar, the electromagnetic emission is not expected to be simply 
described by eq. (4-1) due to the change in geometric configuration and the moment of inertia 
of the pulsar. We now consider the evolution of the electromagnetic emission during spin-up. In 
accordance with our assumption of the flux-freezing, we consider the pulsar magnetic field as a 
dipole field described by 

B(R) = 5£(|_ii^_|l. ,4-33) 

The total flux passing through a closed surface enclosing the pulsar vanishes (i.e. V • B = 0). For 

a chosen rotation axis, 3-axis, wc assume that in the comoving frame the magnetic moment of the 
pulsar lies on the plane defined by the 3-axis and the 1-axis and hence /i^ = (/i*,i, 0, /v,*,3). The 
conservation of the magnetic flux through the 23-plane and V • B = (or the vanishing total flux 
for the closed surface around the pulsar) combine to require that vary as (cf. Finn Sz. Shapiro 
1990) 

R I - fl^Y^^ 1 

'-''^ Ewi-(ls/w =''[i) Ewi-ik/i^r) ^'-''^ 

where i?* = R/ihh)^^^ and E{x) is the complete elliptic integral of second kind (Abramowicz & 
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Stegun 1964). Since the electromagnetic emission from the dipole is given by 



Lem oc 



(4-35) 



we can evaluate how the electromagnetic power varies as the pulsar spins up. Roughly, we expect 
Lem oc Rl^t assuming the magnetic field is a fossil subject to flux- freezing. Thus, the dipole 
power behavior depends on whether or not the angular velocity increases faster or slower than R^. 

Fig. 3(b) shows the resulting exact behavior of the electromagnetic power during the evolution 
driven by the gravitational radiation emission. The electromagnetic power increases roughly by a 
factor ~ 2 which is due to the spin-up and the change of the magnetic moment. That is, during 
the period of the gravitational radiation emission, despite the loss of angular momentum, the 
spin- up of the MSP occurs and Lem remains essentially constant. Therefore, even in this stage, 
the short time scale tcR is not relevant for the characteristic time scale of the GRB durations. It 
is the time scale tEM, which is generally much longer than tcB., which determines the intrinsic 
emission time scale. If real MSPs differ seriously from incompressible homogeneous ellipsoids, our 
results would not be applicable. It would then be difficult to assess the likelihood of spin-up versus 
spin-down during the period of gravitational radiation, and thus whether or not Lem increases 
or decreases in this period. Another complicated issue is the dynamical instability of the rapidly 
spinning neutron stars right after AIC. The new born MSPs could be dynamically stable and 
secularly unstable only if the rotation rates before AIC lie in a narrow region (Chandrasekhar 1969, 
Shapiro & Teukolsky 1983). For a wide range of pre- AIC parameter space, the new born stars 
are likely to be dynamically unstable (Durisen et al. 1986, Williams Sz Tohlin 1988, Houser et al. 
1994). Such dynamically unstable configurations would lose mass and angular momentum on the 
dynamical time scale. Expansion of the mass shed and shock dissipation would result in a nearly 
axisymmetric equilibrium state which is secularly unstable (either due to gravitational radiation 
or fluid viscosity) and evolves into the non-axisymmetric configuration on the dissipation time 
scale. More sophisticated simulations will ultimately be useful in sorting out these complications. 

Given the lack of any alternative realistic calculations, we apply the ideal homogeneous 
ellipsoid MSP model to the discussions on GRBs. Based on the result that the electromagnetic 
luminosity does not change significantly on a time scale ^ tan, we arrive at the following picture 
for the rapidly spinning pulsars from AIC. Initially if the AlC-formed pulsar rotates faster than 
the critical rate, the angular momentum is lost on the gravitational time scale 



and the pulsar spins up. During this stage, newly fromed MSPs could be identified as progenitors 
of GRBs through detection of gravitational radiation. As the equatorial eccentricity ei2 0, the 
evolution of the electromagnetic power occurs on the typical electromagnetic time scale 




(4-36) 




(4-37) 
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For pulsars with initial spin rates below the critical rate, the initial phase is missing and the 
overall electromagnetic emission time scale is given by Iem- As pointed out in Blackman et al. 
(1996) (see below), the observed short duration of GRBs and the large apparent luminosities 
would then require that the emission is sharply beamed in a jet which precesses past the line of 
sight to account for the short time scales. 

The two classes of MSPs (distinguished by whether the spin rates at birth are greater or less 
than the critical rate) could be observationally distinguished by the detection of the gravitational 
radiation. In the faster class, the initial phase would be characterized by the emission of the 
gravitational radiation. This is absent in the slower class. 



4.3. Precession of Gamma-ray Jets 



The overall energy budget requires relativistic beaming (eq. (4-5)). If their energy source is 
electromagnetic radiation, the beamed GRB jets may sustain their luminosities for time scales 
longer than the observed burst durations. This is not inconsistent with observations because these 
jets will precess due to gravitational interaction with their binary partners. Jet precession is a 
second important feature of the present cosmological GRB scenario as it accounts for the fact that 
the observed durations of the bursts are short, even though the jet lifetime in the precessing frame 
can be long (cf. eq. (4-37)). 

Blackman et al. (1996) have pointed out the three important precession frequencies in 
the pulsar binary systems (Thorne et al. 1986): The first is the Newtonian tidal torque 
(J7y). The second is the interaction between the pulsar spin and the gravitomagnetic field 
associated with the secondary star's spin The third (J^g) comes from the pulsar spin 

interaction with the gravitomagnetic field associated with the secondary's orbital motion and 
the secondary's gravitational field, the space-curvature precession, and the spin-orbit precession 
from the gravitomagnetic field induced by the orbital motion of the pulsar in the secondary star's 
gravitational field. The three frequencies, for a maximally spinning neutron star and companion 
are (Blackman et al. 1996, Thorne et al. 1986) 
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(4-39) 



(4-40) 



where fifej^ is the binary orbital frequency. For a sharply beamed jet, the first precession time 
scale above is sufficiently short enough to move the jet out of the observer's line of sight on a very 
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short precession time scale, while the third precession mode would keep the jet from returning to 
the line of sight on > lOOyr time scales. The observed durations (Meegan et al. 1994) could be 
short (~ Is) even when the intrinsic emission time scale 0{W)yr is much longer and the absence 
of repeaters would be accounted for by the third precession mode (Blackman ct al. 1996). 

In the present picture, the source of emission is the electromagnetic dipole emission. For 
a wide range of parameters, the precession would be required if the emission is beamed. The 
relativistic beaming is naturally expected in the cosmological scenario. Even when the precession 
is not required due to intrinsically short emission lifetime (cf. eq. (1-8)), the precession is a 
natural outcome of binary accretion systems. The precession always tends to give apparent short 
durations and to ensure absence of repeaters. 

The apparent duration is determined by the first precession mode (eq. (4-37)) 

rrf,,~(rfiG)"\ (4-41) 

where T is the bulk Lorentz factor for outflowing gamma-ray emitting blobs. The repetition time 
is constrained to satisfy (eq. (4-40)) 

Trep ~ 27rO^^ > lo'^sec, (4-42) 

because of the absence of repeaters. We can combine these two relations with (4-38) and (4-40) to 
obtain a lower limit on the binary radius and a lower limit on F. Using the above parameters for 
the mass scalings, these limits are given by 

D > 7.3 X 10^ cm (4-43) 

and 

F™„ > 4 X 105(Trf„,/Ls)-i. (4-44) 

In the next section we describe a process which gives such a Lorentz factor. (For SS433, Martin & 
Rees, 1979 used the observation of 1 precession period to constrain the binary companion star.) 



4.4. Synchrotron-Inverse Compton Mechanism for Gamma-Ray Emission 

The details of the gamma-ray emission mechanism still remains an outstanding issue. As 
a direct application of the derived spin-field correlation, we consider a particular gamma-ray 
emission mechanism outlined in Blackman et al. (1996). The derived spin-field correlation for the 
AIC pulsars limits the available field strength for a given spin rate, which could set an interesting 
constraint on gamma-ray emission mechanisms relevant for MSPs. 

We consider the synchrotron-inverse Compton mechanism (Asseo et al. 1978). The charged 
particles (electron-positron pair plasma) are accelerated to large F's (which is essential to avoid 
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the run-away pair production [Krolik k. Pier 1991, Yi 1993] that would make the emitting material 
optically thick) by large amplitude electromagnetic waves propagating outside the radius (Usov 
1994) 

i^//-[10-cH(^) (^) (4-45) 

where the flux-freezing and force-free conditions are broken as the charge density decreases below 
that which can sustain the Goldreich-Julian density (Usov 1994, Goldreich & Julian 1969). At 
R = Rff, the number density of the pair plasma becomes 



The electron acceleration parameter cr// at this radius becomes 

where F// = AirR'jfCNff is the number flux at i? = Rff (Michel 1984, Usov 1994). Solving the 

equations of motion for a particle in a pulsar wind zone subject to electromagnetic forces, leads to 
the result that relativistic electromagnetic waves of frequency fi* can accelerate pair plasma to 

i^max ~ (^ff ~ 10 — TT— — TT^T— T (4-48) 



KlO^cmJ yW^Gj \10^s-\ 

(Asseo et al. 1978, Michel 1984). As mentioned earlier, the emission is beamed into an angle 
~ Tmax correspondingly large luminosities received per solid angle may be the result of 

this relativistic beaming (Yi 1993, Usov 1994, Blackman et al. 1996). 

As in synchrotron radiation (e.g. Rybicki and Lightman, 1979), the characteristic emitted 
frequency of the optically thin synchro-Compton radiation is proportional to r^^^, or a'jf from 
(4-48). In particular (Asseo et al. 1978) 

The synchrotron-Compton scattering spectrum's tail can extend up to 



4/3, 



.Bff/^.4 - [10-.-] {^f {^f" {-^f (4-50) 



for a dipole-like field which at R = Rff has B{R = Rff) = Bff. 

Observable gamma-ray emission from the synchrotron-inverse Compton scattering mechanism 
requires 

\ / R*: \ ( B^ 



> 40 (4-51) 
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or for = 10^ cm 



which can be compared with the spin-field correlations derived earlier. 

Fig. 4. schematically shows the two spin-field correlation lines for the two magnetized disk 
models and the gamma-ray emission line derived above. Phenomenologically, the two models 
differ only in the effective strength of the field in the accretion disk S. That is, for cc = 0.1 and 
7 = 1 or 5 ~ 3 and for ^ymax = 1 or 5 = 1 and the difference between the two models is mainly 
due to the effective field strength S. Given the fact that 5 > 1 is most likely, the condition for the 
gamma-ray emission results in a quite stringent constraint on the spin-field values. It is very likely 
that gamma-ray emitting pulsars from AIC will have < a, few xlO^^G and 17* > lO^s^^ which 
puts most of the gamma-ray MSPs in a regime close to that of the fizzlers. However, taking into 
account the uncertainties in S, it is likely that there are a small number of gamma-ray MSPs with 
~ a few xlO^^G and f2* ~ 10^s~^ with gamma-ray emission. Any realistic pulsar gamma-ray 
emission mechanisms are expected to be constrained in a similar fashion. 



4.5. Very Strongly Magnetized Pulsars and Observed Gamma-Ray Bursts 

If the gamma-ray emission is isotropic and burst sources are cosmological, > 10^^ G dipole 
fields are required (eqs. (1-3), (1-7)). Since the flux-frozen fossil field is limited below ~ lO^^G in 
MSPs formed by AIC, strong fields > 10^^ G would have to be generated in situ. 

A traditional a — J7 dynamo can in principle generate mean fields exponentially through a 
combination of helical turbulence, diflFerential rotation, and turbulent diffusion (e.g. Parker, 1979; 
MoflFatt, 1978). Duncan &; Thomson (1992) suggest that for young MSPs, the available turbulent 
energy and shear motions, which allow wrapping of field lines around the star, can generate dipole 
fields of order 10^^ G. 

We will see below why this field strength may still be too small, but even so, there are other 
challenges to invoking in situ field generation in MSPs: In the Galaxy for example, exponential 
dynamo growth of magnetic field is essential because of the long growth time (e.g. Ruzmaikin 
et al., 1988) relative to the total available Galactic lifetime. However in the case of MSPs, even 
if the full feedback leading to exponential mean field dynamo growth were inoperative and only 
linear growth by differential rotation (e.g. Ruzmaikin et al., 1988) were present, the short rotation 
periods would imply that strong dipole fields should still be generated if the shear energy were 
available to stretch the field lines by flux freezing. 

However, it is precisely the ease with which this could be expected to occur within the 
assumption of flux freezing that might pose a challenge. If the differential rotation energy of 
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MSPs (or sub-millisecond pulasars) could be converted into magnetic fields (Duncan & Thompson, 
1992), then most MSPs would have ~ 10^^ Gauss fields, or no MSPs should be observed: If MSPs 
are formed at birth, then these objects should spin down rapidly (~ Ihr) and would no longer be 
seen as MSPs. If the objects are formed later by spinning up, then the dynamo should generate 
strong fields at this stage, and and here again the objects should spin down. One difference is that 
the latter type of objects would have a crust, which may somewhat keep the field from escaping as 
a dipole, but the extent to which this would keep the rotation from spinning down is not entirely 
clear. 

If MSPs are formed by AIC then the absence of such strong fields suggests that such field 
growth is not a canonical feature of MSPs. If there were such AIC formed MSPs, then we would 
expect to see Galactic sources with available dipole power of > lO^'' ergs/sec without supernova 
precursors. Observed field strengths of MSPs are much lower than the dynamo field strengths of 
~ IQ-'^^G of Duncan &: Thompson (1992). In any case, complications with invoking in situ field 
generation suggest that constraints on spin-magnetic field relations from just AIC considerations 
may be the most important in constraining MSP properties. 

Let us suppose, however, that strong dipole fields > lO^^G could be formed in MSPs, (for 
example by a mechanism that induces dynamo growth only in newly formed pulsars) and address 
their consequences in the cosmological scenario. A MSP with 17* ~ 10'^s~^ would provide a large 
enough isotropic luminosity 2 x 10^^ erg/ s only when the field strength 

> [10^6G](^/0.1)-V2 (4.53) 

which depends on the gamma-ray efficiency (e.g. Usov 1994). Even when 
Lem ~ [10^^er5/s](^/0.1)~^, the expected duration of the event is estimated at Iem ~ 5s 
which would account only for the longer burst durations only. This constraint is based on 
our result from section 4.2 that spin-up (rather than spin-down) occurs during the period of 
gravitational radiation, and that Lem therefore does not decrease on a time scale t^/j <^ Is. More 
realistic extensions of the calculations in section 4.2 are of great interest but they are beyond 
the scope of the present investigation. This implies that short bursts with r^ur ~ 0.1s would 
require ~ lO^^G dipole fields, which is highly unlikely. Thus in order to account for both 
cosmological luminosities and short durations, the field strength has to be close to i?* > lO^^G. 
On the other hand, if gamma-ray emission is relativistically beamed with F > 10 as often seen in 
extragalactic jets, ~ 10^^(r/10)~"'^G would be sufficient to account for the observed luminosity. 
In this case, the overall duration would become as long as ~ 50s (eq. (1-8)) which is certainly 
too long to account for short duration bursts. If F S> 10, even though precession would make 
the apparent duration short, the observed luminosity would be too high unless the gamma-ray 
efficiency ^ ~ F~^ is very low. Such a low electromagnetic efficiency (and hence gamma-ray 
efficiency) is not surprising in GRB models based on MSPs as most of pulsar rotational energy is 
carried away by gravitational radiation. 
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5. Summary and Discussions 

We have examined the formation of MSPs from electromagneticahy quiet AIC of magnetized 
white dwarfs. Using magnetized accretion disk models, we have derived a relation between the 
stellar spin and magnetic field, the latter of which is assumed to be flux-frozen during accretion 
and collapse. For the high mass accretion rates required for AIC, the MSPs are born with a typical 
field strength less than ~ 10^^ G whether or not they reach their equilibrium spin state during the 
pre-coUapse stage. A significantly different spin-field relation is expected only when the secondary 
star is a dwarf star with a very small binary separation, and then the spin and binary orbits would 
be synchronized. Spin-binary orbit synchronization would be absolutely required for formation of 
the strongly mag nctizcd (~ 10^^ G) MSPs but the derived spin-magnetic field correlation strongly 
suggests that rapid spin and very strong magnetic fields are mutually exclusive. 

This conclusion remains valid within the range of initial conditions we adopted, which 
include an almost entire range of initial white dwarf spin periods. This is not surprising given 
the short spin-up time scale (eq. 3-13) leading to near spin equilibrium. Siginifcantly different 
initial conditions (such as white dwarf initial spin and mass) and physical conditions (white dwarf 
magnetic field, mass accretion rate), which are still appropriate for AIC, in principle, could affect 
the spin-field correlation. Within our disk-magnetosphere model, however, it is highly unlikely 
to have MSPs with ^ lO^^G. The spin-orbit synchronization however remains a possible 
route to MSPs with very strong mag netic fields > lO^^G (Usov 1992). It wih be crucial for the 
synchronization model to detect ~ lO^G white dwarf field. Alternatively, it is conceivable that 
the disk-magnetosphere interaction is substantially different from what we have described. In 
the presently available alternative models, it has been argued that the spin-up by magnetized 
accretion disks is much less efficient than in the present model (e.g. Shu et al. 1994, Lovelace et 
al. 1995, Li 1996). In this sense, our spin-field correlation can be regarded as most favorable to 
formation of MSPs with strong flux-frozen magnetic fields and hence as providing an upper limit 
(~ 10^2 Q) foj. the MSP frozen field. 

Given the required high mass accretion rates for AIC, and the absence of the observed 
~ lO^G field white dwarfs, any ultra-strong dipole magnetic fields ~ lO-'^^G MSP are extremely 
difficult to produce by AIC (cf. Usov 1992) and, as discussed in section 4., they do not seem to 
be characteristic features of Galactic MSPs as would be expected if the in situ field production 
(Duncan & Thompson, 1992) were a canonical feature of MSPs. A small fraction of magnetized 
white dwarfs could have ~ 10^ G although such systems have not been observed. These systems 
would become MSPs (via AIC) only when the spin-orbit synchronization occurs prior to AIC. 
This is because the disk-magnetosphere interaction would most likely to give long spin periods 
(eq. 3-7). Even if such MSPs existed, the constraints from the observed luminosity and duration 
severely limit the available parameter space which essentially requires 5* > 10^^ G for bursts with 
shortest durations. 

The spin-field correlation puts an upper limit on the intrinsic electromagnetic power which is 
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fax less than the observed luminosities of the cosmological GRBs if the emission were isotropic. We 
conclude that the gamma-ray emission has to be strongly beamed in an AIC scenario. The small 
required electromagnetic power would lead to GRBs much longer than observed if it were not for 
precession of the beamed gamma-ray emitting jets. The typical binary precession periods could 
account for the observed durations (Blackman et al. 1996) and absence of repeaters. Anisotropic 
emission, mediated by magnetic fields in an AIC scenario is consistent with overcoming the baryon 
contamination problem thoroughly investigated by Ruffert et al. (1996). 

Adopting the synchro-inverse Compton mechanism for the gamma-ray emission as an example, 
we have applied the condition for detectability of the gamma-ray emission from MSPs. This 
mechanism favors gamma-ray emission from pulsars with relatively weak fields 10^ ^G) but fast 
spins. These pulsars are likely to have initial conditions similar to those of the pre-collapse white 
dwarf systems which result in the fizzlers. In this sense, the fate of the fizzlers is quite interesting 
(Tohlin 1984) as is the exact fate of the dynamically unstable new born stars from AIC (e.g. 
Houser et al. 1994). 

For rapidly rotating pulsars, whose spin frequencies are larger than the critical frequency, the 
initial phase of the emission is dominated by the gravitational radiation. In the simplified picture 
of the homogeneous ellipsoids, during the initial phase, the electromagnetic emission remains 
nearly constant or increases despite the loss of the angular momentum. After the short initial 
phase < Is, the pulsars become axisymmetric, the gravitational radiation becomes small, and the 
emission is dominated by the electromagnetic emission. For initially slowly rotating pulsars, the 
initial gravitational radiation dominated phase is absent. In both cases, the binary precession 
is crucial for the short durations. Hence, apart from the initial phase, the long electromagnetic 
emission phase is the inevitable outcome for the pulsars formed in AIC. The existence of the initial 
phase in the MSPs spinning above the critical rotation frequency could be distinguished by the 
accompanying gravitational radiation. 

Blackman et al. (1996) estimated the required AIC event frequency for the GRBs. The 
observed burst rate ~ 10~^yr~^ per galaxy requires the source event rate ~ per galaxy. 

Though the estimated rate is indeed larger than observed rate because of the beaming, the 
discrepancy is not nearly as extreme as it would be if the jet lasted only a duration time: The 
long lifetime of the emission from any given object means that the probability for observing that 
object is not merely the beam angle but the angle swept by the beam that can pass through the 
line of sight during the lifetime of the object. The production rate combined with the long life 
time of the MSPs suggest that there would be < 10 jetted objects in a typical galaxy. When the 
jet beams are not oriented toward us, these objects would appear similar to the normal pulsars. 
The known ~ 10^ cataclysmic variables (as progenitors of AIC) in our Galaxy and the ~ 10^ 
low mass X-ray binaries (related to MSPs) then imply that a fraction of AIC events may have 
resulted in gamma-ray emitting MSPs with the characteristic intrinsic luminosity ~ 10^^~^^ erg'/s. 
It is interesting to point out that the expected events from electro-magnetically quiet AICs with 
the canonical luminosity ~ W^^~'^'^erg/s (cf. eq. (4-4)) and the time scale ~ 1 — lOyr (cf. eq. 
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(4-6)) from new-born MSPs have not been identified. In our scenario, the observed GRBs (with 

relativistic beaming and precession) are the results of such AIC events. We have shown that these 
objects are not the canonical result of AIC; they fit into a small but available parameter space, 
and would not ovcrpopulatc the MSP population. 

Only if the dipolc magnetic fields were > 10"*^^ G would the luminosity of the pulsar be high 
enough to make an isotropic ORB emission at a cosmological distance and these fields would have 
to be produced by neutron star dynamos, if these are opertative. 

If the gamma-ray emission region is close to the neutron star, the critical field strength 
> 4.4 X 10^3 G in the gamma-ray emission region could produce observable features in gamma-ray 
spectra through photon splitting (e.g. Adler 1971, Baier et al. 1996, Adler & Schubert 1996) or 
photon-magnetic field attenuation (e.g. Ho et al. 1990). If the emission region is far away from 
the neutron star surface, however, such spectral signatures are not expected. If the GRBs are 
cosmological, the severe requirements due to the gamma-ray transparency (Krolik & Pier 1991) 
and baryon loading problem (Meszaros &: Rees 1996) point to the gamma-ray emission region 
well away from the neutron star surface. The field strength in the gamma-ray emission region is 
likely to be lower than the critical field strength, which makes the direct confirmation of the field 
strength > IQ-'^^G unlikely. 

We thank G. Field, P. Kumar, R. Narayan, M. Rees, and E. Vishniac for many helpful 
suggestions, discussions and comments during various stages of this project and the referee of our 
previous paper who raised the issue of the electromagnetic emission from rapidly spinning pulsars. 
Helpful suggestions and criticisms by the referee, V. Usov, are gratefully acknowledged. I. Y. 
acknowledges support from SUAM Foundation. 
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Fig. 1. — Torque exerted on the star by the magnetized disk as a function of the disk truncation 
radius. The upper curve corresponds to the model where the pitch of the azimuthal component 
of the magnetic field is limited by the nearly force-free condition in the stellar magnetosphere. 
The lower curve corresponds to the model with the diffusive field loss mechanism. The difference 
is mainly due to the torque contribution from the region outside the corotation radius where the 
two models differ. The equilibrium spin is achieved at x^g = 0.9148 in the lower curve and at 
Xeq = 0.9502 in the upper curve. 
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Fig. 2. — Spin — magnetic field correlation after collapse from white dwarfs to neutron stars. Initial 
mass (Mj) and mass accretion rate are shown in the panels. Initially all white dwarfs are taken as 
sufficiently slowly rotating with varying magnetic field strengths. The thin lines correspond to the 
pitch-limited magnetic torque and the thick lines correspond to the model where the field loss is 
due to turbulent diffusion. In each panel, the dashed line corresponds to the hypothetical critical 
rotators which are assumed to collapse with the conditions exactly at critical rotation. The dotted 
line corresponds to the fizzlers and the solid line to those directly collapsing to neutron stars. 
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Fig. 3. — Time evolution of a rapidly spinning neutron star approximated as a Jacobi ellipsoid. 
The panel (a) shows the spin-up of the pulsar as a result of the gravitational radiation emission. 
The panel (c) shows the evolution of the equatorial eccentricity, ei2, defined in the text. The star 
evolves toward an axisymmetric spheroid with a non-zero 613. The variation of the gravitational 
radiation luminosity (normalized by the luminosity at r = 0) is shown in panel (d). The panel (b) 
is the evolution of the electromagnetic luminosity (normalized by the luminosity at r = 0) as a 
result of the gradual change in the moment of inertia and loss of angular momentum through the 
gravitational radiation. 
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Fig. 4. — Schematic constraints for gamma-ray emission from pulsars. The gamma-ray emission 
from the synchrotron-inverse Compton scattering mechanism is possible when the parameters are 
above the solid line i.e. the spin has a lower bound defined by the solid line for a given field 
strength. The two dotted lines are the predicted spin-magnetic field correlations marked by the S 
values. The upper dotted line is for the pitch-limited model with S = 1 and the lower dotted line 
is for the diffusive loss model with 6 = 3. 



